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Abstract 



We consider a holographic model of QCD, obtained by a very simple modification 
of the original construction, which describes at the same time the pion mass, the QCD 
anomalies and the baryons as topological solitons. We study in detail its phenomeno- 
logical implications in both the mesonic and baryonic sectors and compare with the 
observations. 



1 Introduction 



The virtues of 5d theories as phenomenological descriptions of the QCD hadrons in the 
large-iV c expansion have been widely discussed in the literature [1-8] and can be briefly 
summarized as follows. From the theoretical point of view, these models resemble large- N c 
QCD in that they contain infinite towers of weakly interacting mesons while the baryons 
are solitons, similar to the Skyrmions of the non-linear a-model. Differently from ordinary 
Skyrmions, however, the 5d soliton size remains finite in the weak coupling (large- N c ) limit 
and is therefore parametrically larger than the length scale at which the 5d model enters 
the strongly coupled regime. The 5d Skyrmions are genuinely macroscopic objects, whose 
physics is perfectly within the reach of the effective theory as it must also be the case for 
the baryons in the true theory of large- N c hadrons [9] . 

For what concerns phenomenology, it is already quite remarkable that the holographic 
implementation of the QCD chiral symmetry and of its breaking automatically leads to 
towers of vector and scalar mesons with the quantum numbers of the observed p, u, a±, f\, 
a 0) /o? vr(1300) and ^(1295), but even more remarkable is that, when extrapolating the model 
to the physically relevant case of N c = 3, an agreement at less than 10% with observations is 
found. This success is non-trivial because the model contains an extremely limited number 
of free parameters that can be adjusted to reproduce many observed meson couplings and 
masses, and is mainly due to a set of phenomenologically successful "sum rule" relations 
among the predictions. The latter originate from the 5d structure of the theory and are 
independent of many details of the model, such as the background 5d metric, whose choice 
is therefore unessential for phenomenology as discussed in [10] . The overall agreement is less 
good in the baryonic sector, but still compatible with the expected 1/N C ~ 30% deviations. 
Furthermore, one recovers in a non-trivial way several features of the large- N c QCD baryons 
such as the scaling of the currents form factors with N c , the fact that the isovector electric 
and magnetic radii diverge in the chiral limit (as in QCD [11]) and other large-distance 
behaviors [12]. In holographic QCD, and in the present article, an important role is played 
by the Chern-Simons (CS) term of the 5d action whose presence is required, with a fixed 
coefficient, by the need of reproducing the QCD global anomalies. The CS is therefore 
the 5d analog of the gauged Wess-Zumino-Witten (WZW) term of the standard 4d chiral 
theory and, like the WZW, it is responsible for the it — > 77 decay and for "naive parity" 
breaking in the Goldstone interactions [13]. The CS actually contains the WZW, in the sense 
that it reduces to it in the low-energy description of the Goldstones obtained from the 5d 
theory, but it also contains other interactions such as the well-measured "anomalous parity" 
couplings g p7T -y, g^wy and g upw . As a result of the 5d structure of the theory, these couplings 
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are predicted and are found to be in good agreement with the observations. In the case of 
exact chiral symmetry (and only in that case, as we will show) the CS term is also crucial 
for the physics of the baryons because it stabilizes the 5d Skyrmion solution and makes its 
size scale like a constant for large N c . 

Many of the above-mentioned results, and in particular the ones related with the CS term 
and with the baryons, have however been only established in the minimal holographic QCD 
model with exact chiral symmetry; the aim of this article is to generalize them to the more 
realistic case of explicit chiral breaking (and therefore non- vanishing mass for the Goldstones) 
and to check if the previously outlined general picture survives. In order to complete this 
program we first of all need a model with explicit breaking, but we cannot simply employ 
the original one of [2-4]. A very simple variation is needed in order to incorporate the 
CS term and the 5d Skyrmions that were not considered in the original literature. Like 
in [2-4], our model will be a 5d U(2)l x U(2)r gauge theory with a 5d scalar field $ in 
the bifundamental of the group, whose presence is necessary to parametrize the breaking of 
the chiral symmetry due to the quark mass term. The only difference with [2-4] is in the 
boundary conditions for the gauge fields at the IR brane, where we break the chiral group 
to vector (as in [5]) instead of preserving it and making the spontaneous chiral symmetry 
breaking arise exclusively from the VEV of $. 1 This choice of boundary conditions would 
be a priori equivalent to the original one from the model-building perspective, and costs no 
more parameters, but is necessary if willing to incorporate the CS term. As we will discuss 
in more details in the following section (see also [14]), the reason is that the CS gauge 
variation is an integral on the boundaries of the 5d space and therefore receives contribution 
from both the UV and the IR brane. While the UV term is welcome because it reproduces 
the QCD global anomalies, the IR one must cancel because the IR group is gauged and a 
non-vanishing variation would have on the model the same effects of a gauge anomaly. The 
symmetry-breaking boundary conditions, which are actually equivalent to gauging only the 
vector subgroup at the IR, enforce this cancellation. On top of this, symmetry-preserving IR 
boundary conditions could make the 5d Skyrmion unstable or even the classical 5d Skyrmion 
solution not to exist at all because the baryonic charge 2 

oc f MO - ws(r)] + / MJ)-ws(r)] , (1) 

J r— >oo J z=z m 

1 Actually, another difference with [2-4] is that the gauge group was taken there to be the chiral SU(2) 
while our model is based on U(2). 

2 For the notation, see [6-8] or the following Section. 
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would not be quantized and therefore not topologically conserved. It is indeed possible to 
change continuously the IR contribution to B (the last term in the equation above) without 
affecting the rest by a local variation of the fields at the IR boundary. This is impossible with 
symmetry-breaking boundary conditions because the IR contribution vanishes identically and 
the baryonic charge is quantized as shown in [6-8]. 

The paper is organized as follow. In Section 2 we first of all describe the model which is, to 
our knowledge, the first one containing at the same time the explicit chiral breaking, the QCD 
anomalies (i.e. the CS term) and the 5d Skyrmion. We also discuss its phenomenological 
implications for the physics of vector and scalar mesons, which because of the change in 
the boundary conditions are a priori different from the ones derived in [2-4]. We actually 
find a quite similar phenomenology, the only remarkable difference being in the mass of the 
7r(1300) pseudoscalar meson that was impossible to fit in the original model [4] while it is 
easily accommodated in our case. We obviously also take into account the processes mediated 
by the CS term that were absent in the original model, and eventually perform a complete 
fit of the meson observables which allows us to fix the parameters and to quantify the level of 
agreement with observations. In Section 3 we study the physics of the baryons, generalizing 
the analysis of [6-8] where the chiral case was considered. In that section we check that the 
existence and calculability of the 5d Skyrmion is maintained in the non-chiral case, and see 
how much the predictions change after the deformations required to incorporate the chiral 
breaking in the holographic model. The latter deformation is not mild because the newly 
introduced 5d scalar does not decouple from the vector mesons, and therefore from the 5d 
Skyrmion solution, even when the chiral symmetry is restored. We also compute the isovector 
form factor radii which were divergent in the chiral case and become now finite, as expected 
in QCD, due to the chiral breaking. Finally, in Section 4, we present our conclusions. 



2 Massive Pions and Scalar Resonances 



As a starting point for the description of the model, let us consider the QCD partition 
function Zg C£ ,[Z,r, M\ in the presence of sources for the left and right global currents and 
for the quark mass operator. The latter reads 



Z QCD [\,r,M] = J VVexp iS Q c D [*] +^ J d 4 x Tr (\,f L + r,f R - Ms - s*M*) 



(2) 



where \& collectively indicates the QCD fundamental fields, Sqcd is the massless QCD action, 
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1 J l r^Ql r are the currents of the U(2)l x U(2)r chiral group and Sij = q 3 L q l R is 



the quark bilinear. Due to anomalies, the partition function changes under local chiral 
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transformations as 3 



Zqcd 



\^\v^\g L Mg R ] = e u Z QCD [1, r, M\ , (3) 



where the vector sources transform as gauge fields {i.e., lj^ = ^[1^ + ^(9^]^ and analogously 
for r M ) and .4. stands for the QCD global anomaly which, as discussed in appendix C, can 
be put in the form 

A= ^I H(<*r,> 1 ) -"*(<**>*)] ( 4 ) 

by the addition of suitable local counterterms. In words, what eqs. (2) and (3) mean is that 
2-flavor (large- N c ) QCD is a theory endowed with an U(2) L x U(2) R anomalous global sym- 
metry and with a non-dynamical spurion field M, associated with the quark mass operator, 
in the bifundamental of the group. This symmetry, with its explicit breaking parameter M, 
must obviously be present in any phenomenological description of hadrons such as the one 
we want to construct. 

Following the holographic method, we incorporate the chiral symmetry and the spurion 
by introducing 5d fields with appropriate quantum numbers (Lm, Rm and $) associated 
respectively to the sources 1 M , r M and M, and identify the latter with the values of these 
fields at the UV-boundary. For the gauge fields the model is exactly as in [7], and all what 
we have to add is the scalar <3> which, according to the holographic prescription, transforms 
as 

$ $(fl£,9H) = g L <f>g R , (5) 

under the local 5d U(2) L x U(2) R group and is subject to the following UV boundary con- 
ditions ^ 

*U uv = f— ) M, (6) 



IR 



with A~ defined in eq. (19). The need for the rescaling in M is a technicality associated 
with our choice of the AdSs geometry, and will be explained in the following. For simplicity, 
and since we are considering the two-flavor case in which the vector symmetry breaking is 
negligible, we take the spurion VEV to be of the form 

M = M q l, 



3 The equation which follows is only valid in the large-iV c limit in which the U(\)a~SU(N c ) 2 anomaly 
can be neglected and the rj' only gets its mass from the explicit chiral symmetry breaking. Given that the 
expansion parameter of the 5d theory is interpreted as 1/N C , it is completely correct to stick to eq. (3) at 
the leading order, though it would be interesting to see what kind of next-to-leading order corrections the 
inclusion of the U(1)a anomaly leads to. The problem of incorporating the rj' mass in holographic QCD has 
already been addressed in [15], though its connection with the anomalies and the 1/N C expansion has not 
been outlined. 
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which respects U(2)v- 

The partition function of the 5d model, which we would like to identify with the QCD 
one in eq. (2), is defined as 

Z [1, r, M] = J VL M VR M V<& exp {i S 5 [L, R,$]}, (7) 

where the dependence of the r.h.s. on the sources arises from the boundary conditions on 
the allowed field configurations. The gauge part of the 5d action is given by a kinetic part 

S 9 = ~ j d A x jP " dza{z)^^[L MN L MN ] + l -L MN L MN + {L i?}} , (8) 

and a Chern-Simons part 

Scs = ^- 2 j dS X ^e MNOp QL M Tr[L NO L PQ \ + ^e MN °^L M L NO L PQ - {L » i?}} , 

(9) 

having parametrized the fields as = L^a a /2+L M l/2 and analogously for R M . As discussed 
before, the UV boundary conditions for the gauge fields are 



J fl | Z — Zuy 



= 1 M , R M \ z=Zvv = r M . (10) 



The action for the scalar is 



S <S , = M 5 j d"x r dza\z) {Tt[(D m ^ D m §\ - a(z) 2 M 2 Bulk Tr [$t$] j , (H) 



Zuv 

where we defined the covariant derivative 



D M <5> = d M <5> - i L M $ + i $R M . (12) 

If a local 4d chiral transformation is performed on the sources, as in eq. (3), the UV boundary 
conditions get modified but this change can be reabsorbed in a change of variable in the path 
integral of the form of a local 5d gauge transformation g^R = exp(iaL,n) which does not 
reduce to the identity at the UV brane. The kinetic part of the gauge action S g and the 
scalar action S$ are invariant under this transformation, while the CS part (as discussed in 
appendix C) gives the variation 

AaScs = 24^ ( fM( a L,l)-w\(a R ,T)] - f[u\(a L ,L)-Uj\(a R ,R)] \ . (13) 

\J Z — Zuy J Z — Zir / 

The UV contribution in the above expression is welcome, because it coincides with the 
anomalous variation in eq. (3). On the other hand, the IR terms in A a Scs have no coun- 
terpart in 4d QCD and must be cancelled. The chiral-breaking conditions 

(I> " R„) = , (L M5 + R M5 ) \z=z m = (14) 
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enforce this cancellation, which, on the contrary, would not occur in the case of the symmetry- 
preserving Neumann conditions considered in [2-4]. The choice (14) of boundary conditions 
is motivated not only by the need of correctly reproducing the QCD anomaly, as the previous 
discussion shows, but also by the consistency of the theory in the presence of the CS term. 
We can introduce a term in the action only if it respects all the symmetries which are 
gauged. As can be seen from es. (13), the CS term in eq. (9) is not invariant under local 
axial transformations which do not vanish at the IR boundary, thus it can be consistently 
introduced only if the local axial symmetry at the IR boundary is not gauged, as implied by 
the conditions in eq. (14). 

For the scalar we impose the IR boundary condition 

$| z=2lR = ±£1, (15) 

where £ is a positive parameter and fields satisfying eq. (15) with both signs, associated 
to two disconnected sectors, coexist in the theory. These boundary conditions (with their 
sign ambiguity) can be thought to originate (as in [2-4]) from an IR localized potential of 
the form V = A ((Re Tr[<£>]) 2 - 4£ 2 ) 2 + rj (2Tr[$t$] - (Re Tr[$]) 2 ) . The latter would enforce 
eq. (15) in the A, 77 — > 00 limit. When dealing with small field fluctuations around the 
vacuum, the subtle sign ambiguity is irrelevant, because only one of the two sectors contains 
a stable stationary point and the dynamics entirely takes place around it. This is because of 
the $ — > — $ symmetry of the lagrangian, which implies that a simultaneous sign change of 
M q and of the boundary condition (15) is unobservable. The pion mass squared is therefore 
m 2 ~ ± M q t; f(zm£)+0(Mq) which, depending on the sign of M q , can be positive or negative 
implying that only one of the two sectors can be tachyon-free. For M q > the vacuum is in 
the "plus-sign" sector and the physics of mesons, described by small fluctuations around the 
vacuum, is totally insensitive to the presence of the other (minus-sign) sector. We will show 
in the following section that the baryons live, on the contrary, in the minus-sign sectors so 
that the sign ambiguity in eq. (15) has to be maintained. 

Having described the model, we can now study its phenomenological implications for 
the physics of mesons, with the aim of comparing it with the observations and of fixing 
its parameters. With respect to the massless case of [7], three new parameters {M q , £ and 
^Buik) have been introduced, but also new predictions can be extracted from the model. 
The scalar $ describes scalars and pseudo-scalars with the isospin quantum numbers of, 
respectively, the a (980), the / (980), the 7r(1300) and the 77(1295); we will use their masses, 
plus of course the pion mass, to fit the new parameters. 4 

4 In our model the axial field Am gives rise to towers of resonances with the quantum numbers of the /1 
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To study the meson spectrum, we have first of all to find the vacuum configuration, which 
is given by the solution of the bulk equation of motion for $ 

D M (a\z)D M <$>) = -a\z)M 2 Bulk <$> , (16) 

for vanishing 4d momenta and with the boundary conditions in eqs. (6) and (15). The result 
is 

m ^j±y\ c (±y S „ W1 , (17) 



with 

( 7 2a 

(18) 



>y 

IR UV 



where we normalized the warp factor as a(z) = z m /z and we defined 



A ± = 2±^/4 + M| = 2±a, (19) 

with M| = z 2 R M Bulk . To obtain a real value for A ± and to avoid a singular behavior of 
the z A piece in eq. (17) at the UV boundary in the z vv — > limit, we must impose the 
constraint —4 < Mj < 0. Notice that it is only because of the rescaling in eq. (6) that we 
obtain a finite VEV in the z VY — > limit while keeping all the other parameters finite. The 
rescaling has allowed us to define a finite quantity, M q , which controls the departure from 
the chiral limit. 

2.1 The Mesons Wavefunctions 

In order to study the properties of the mesonic sector, it is useful to rewrite the scalar field 

as 

$ = (vl + 5) e iP/v , (20) 

where S and P are Hermitian matrices of respectively scalar and pseudoscalar fields. Under 
the unbroken U (2)y symmetry S and P can both be decomposed as 1+3. Form the boundary 
conditions on the $ field in eqs. (6) and (15) we immediately find the boundary conditions 
for the S and the P fields 

^ I Zuv ^ I Zir ^ I Zxjy ^ I Zir ^ ( ^ 1 ) 

and of the rj mesons. These states are predicted to be degenerate in mass with the corresponding mesons of 
the ai and the tt tower. Although the masses predicted for the /i (1285) and for the 77(1295) are close to the 
experimental ones, we decided not to include these predictions in the fit because, in principle, they could 
receive sizable corrections when a proper treatment of the U{1)a anomaly is included (see footnote 3). 
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A convenient gauge fixing choice for studying the properties of the mesons is the gauge 
which eliminates the mixing between the gauge bosons A^, and the scalars A 5 , P. This 
gauge is obtained by introducing in the Lagrangian the terms 

(22) 



£V = 2M 5 a( Z ) Tr 



'9f £ 



V 



a{z) 



d^A" - ^d z (a(z)A 5 ) - Ua 2 (z)v(z)P 



(23) 



where Vm = (Lm + R-m)/2 and Am = (Lm — Rm)/2- The wavefunctions and the masses of 
the mesons can be determined by solving the quadratic equations of motion for the 5d fields 
and imposing the appropriate boundary conditions. The computation is analogous to the one 
described in [3,4], so we will skip here most of the details. The physical degrees of freedom 
are easily identified using the unitary gauge, which corresponds to the limit ^a,v — > °o. In 
this gauge the fields A 5 , V5 and P satisfy the constraints 

d z (a(z)V 5 ) = , P = -^—d z (a(z)A 5 ), (24) 

and the pion field can be identified with the lightest KK mode of A 5 . 

The vector mesons, namely the p and the u mesons and their resonances, are described 
by the KK modes of the vector gauge field V^. In the unitary gauge the V 5 component of 
the vector gauge field is forced to vanish due to the boundary condition V 5 | 2ir Zm = 0, while 
the 4d components satisfy the bulk equation of motion 

d v W» ~ d^V v - -±-d z (a(z)d z V„) = , (25) 
a{z) 

and the boundary conditions 

VJ =0, d z V,\ =0. (26) 

As expected, the equation of motion for the vector mesons is independent of the VEV of the 
5d scalar $, and the masses and the wavefunctions for these states are exactly the same as 
in [7]. 

The equation of motion for the scalar field S is independent of v as well 

d^S - ^d z {a 3 (z)d z S) + a 2 (z)M 2 Bulk S = , (27) 

and the corresponding boundary conditions are given in eq. (21). The KK modes of the S 
field are interpreted as the parity-even scalar mesons, whose lightest states are the ao(980) 
and the Jo (980). 
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On the other hand the axial fields, which are interpreted as the ai(1260) and /i(1285) 
mesons and their heavier resonances, are sensitive to the symmetry breaking induced by the 
VEV of $, as can be seen from the equation of motion 

dv&Ap - dpVAv - -^d z (a(z)d z A^) + 2a 2 {z)v 2 {z)A il = . (28) 
a{z) 

The boundary conditions are 

A4 uv = A4 IR = o. (29) 

Finally, the pion and the other pseudoscalar mesons are described, for finite £4, by 
the fields A 5 and P. These fields mix in the quadratic Lagrangian and it is not totally 
straightforward to obtain their equations of motion in the unitary gauge limit £4 — >■ 00. If 
one simply removes the P field from the Lagrangian through the constraint in eq. (24) and 
then varies the action with respect to A 5 , one obtains indeed the bulk equation of motion 

V [d^A b + 2a 2 (z)v 2 (z)VA 5 ] = , (30) 

which is of fourth order in the z derivatives, having defined the differential operator 

V=l-d z ( - „. \ r d g a(z)) . (31) 



2a 3 (z)v 2 (z) 

From the boundary conditions on P and from the gauge fixing constraint one also gets the 
boundary conditions 

9,(0(^)1^ = 9,(0(^)1^ = 0, (32) 

which however are only two while four conditions would be needed for the fourth order bulk 
equation (30) to have a unique solution. The two remaining boundary conditions cannot be 
obtained directly in the unitary gauge, they arise from the equations of motion at finite £a 
by carefully taking the £a — > 00 limit. From this procedure we recover the bulk equation 
(30) and the boundary conditions in (32), plus the additional constraints 

d z {a{z) [d^ + 2a 2 {z)v 2 {z)V] A,}^ = . (33) 

By the above equation one can prove that the bulk equation (30) can be replaced by the 
reduced second-order one 

d^A 5 + 2a 2 {z)v 2 {z)VA 5 = , (34) 

with the boundary conditions given in eq. (32). The pseudoscalar spectrum, and in particular 
the pion, its lighter state, obviously depends on v(z) and therefore on M q and £. We have 
checked that for positive M q (and £) all the squared masses are strictly positive, while for 
M q < the lightest state is tachyonic as anticipated in the discussion below eq. (15). 
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Figure 1: Mass of the 7r(1300) meson (left) and ai(1260) meson (right) as a function of £L for 
L _1 = 320 MeV and l?M\ ulk ~ —3.8. The shaded band shows the experimental values m^-moo) = 
(1300 ± 100) MeV and m ai(1260 ) = (1230 ± 40) MeV. 

2.2 Fit on the Meson Observables 

Before performing the complete fit, an estimate of the five parameters of our model (L = 
z m — Zw, M 5 , M BuU ,, £ and M q ) can be obtained as follows. The length of the compact 
dimension is fixed by the p meson mass through the relation m p ~ 2.4/ L, from which 
we derive L~ x ~ 320 MeV [3]. The bulk mass parameter is determined by the mass of 
the ao (980) meson. The relation between these two quantities can be derived by solving 
analytically the equation of motion for S (eq. (27)) and then imposing the corresponding 
boundary conditions. This leads to an approximate relation 

A+ 

m ao (980) - 1 - 3 ^ ' ( 35 ) 

which gives A + ~ 2.4 or, equivalently, L 2 M\ ulk ~ —3.8. An estimate of the value of the £ 
parameter can be obtained from the values of the masses of the 7r(1300) and of the ai(1260) 
mesons, which we determine numerically using the previously determined values for L and 
Msuik- As can be seen from the plots in fig. 1, fitting the value of the ai(1260) mass points 
towards the region of small £, £L < 2, while reconstructing the 7r(1300) mass favours larger 
value of the parameter 2 < £L < 3. Combining the two regions we get as a rough estimate 
£L ~ 2. The values of M5 and M q can be extracted respectively from the pion decay 
constant f n and the pion mass m n . These two quantities can be determined by computing 
the holographic Lagrangian for the pion field as explained in appendix A. We find that the 
series expansion for small £L is given by 

fl - ^ (l + A+( ( ltr_ 1} ) > rni - 2(A+ - 2)£M, . (36) 
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135 MeV 


134 MeV 


0.6% 




135 MeV 


133 MeV 


1.2% 


m 7r(1300) 


1300 MeV 


1230 MeV 


5.6% 


m 7r(1300) 


1300 MeV 


1207 MeV 


7.7% 


m p 


775 MeV 


783 MeV 


1.0% 


m p 


775 MeV 


842 MeV 


8.6% 




782 MeV 


783 MeV 


0.1% 


m w 


782 MeV 


842 MeV 


7.7% 


m ai(1260) 


1230 MeV 


1320 MeV 


7.6% 


m ai(1260) 


1230 MeV 


1387 MeV 


12.7% 


m a (980) 


980 MeV 


1040 MeV 


6.5% 


m a (980) 


980 MeV 


1099 MeV 


12.1% 


m /o(980) 


980 MeV 


1040 MeV 


6.5% 


"Vo(980) 


980 MeV 


1099 MeV 


12.1% 


u 


92 MeV 


89 MeV 


3.6% 


u 


92 MeV 


89 MeV 


3.1% 


fp 


153 MeV 


149 MeV 


2.7% 


fp 


153 MeV 


158 MeV 


3.1% 


fu, 


140 MeV 


149 MeV 


6.4% 


f u 


140 MeV 


158 MeV 


12.7% 




6.0 


4.89 


22.7% 


9pnn 


6.0 


5.33 


12.5% 




0.72 


0.71 


1.1% 


9uiTTf 


0.72 


0.74 


3.2% 


Qpirj 


0.22 


0.24 


7.9% 


9fmi 


0.22 


0.25 


12.6% 


9u>pTT 


15.0 


15.6 


3.7% 


9uipw 


15.0 


16.8 


11.8% 


RMSE 






7.7% 


RMSE 






9.6% 



Table 1: Meson observables used for the fit of the microscopic parameters. The table on the left 
shows the results of a fit which minimizes the overall root mean square error (shown in the last 
line). The values of the parameters obtained in this way are L _1 = 325 MeV, M5L = 0.014, 
L 2 M 2 Bulk = -3.7, £L = 2.1 and M q = 31 MeV. The table on the right gives the list of results 
obtained by minimizing the largest deviation from the experiments. With this procedure we get 
L- 1 = 350 MeV, M 5 L = 0.014, L 2 M\ ulk = -3.8, £L = 1.5 and M q = 40 MeV. 

Using the experimental values 5 fl = 92 MeV and m n = 135 MeV, we get M 5 L ~ 0.015 and 
M q ~ 35 MeV. 

A more precise determination of the microscopic parameters can be obtained by per- 
forming a fit on a larger set of well measured mesonic observables. This procedure provides 
also a way to estimate the level of agreement of the model with the experimental results. 
As a simple fitting procedure, we chose to minimize the root mean square error (RMSE) 
of our predictions with respect to the experimental data (see [10] for details on the fitting 
procedure). We remark that for our analysis we take into account only the deviation of 
the theoretical predictions of our model from the central value of the experimental results. 
A more refined procedure, which however would be beyond the scope of our work, should 
also take into account the experimental error with which the various observables have been 
measured. 

The list of observables used in the fit includes the masses of the lightest mesonic reso- 

5 We did not include the electroweak correction in the determination of the observables. For this reason we 
compare the predictions of our model with the experimental value for pion decay constant with subtracted 
electroweak contributions [16] and with the mass of the n°, whose electroweak corrections are negligible [17]. 
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nances as well as some of their decay constants and couplings. The predictions of the model 
and the experimental values are shown in the list on the left of table I. We found that 
the best agreement with the data is obtained for the following values of the parameters 
L- 1 = 325 MeV, M 5 L = 0.014, L 2 M 2 Bulk = -3.7, £L = 2.1 and M q = 31 MeV, which are 
close to the previously estimated ones. The overall agreement with the experimental data 
is quite good and almost all the observables show a deviation from the experimental values 
smaller than 8%, resulting in a RMSE of 7.7%. A somewhat surprising result of the fit is 
the fact that only one observable, namely the g pn7T coupling seems to have a large deviation 
from the experiments. This deviation is however not completely unexpected. In the class of 
models we are considering one gets an approximate tree-level relation [3] 

m 2 p * 3f 2 g% w > ( 37 ) 

which differ by a factor 2/3 form the experimentally well satisfied KSRF relation m 2 ~ 
2f 2 g 2 fmir In our fit the predictions for the pion decay constant and for the p meson mass are 
in excellent agreement with the data, thus the g pn7T coupling must necessarily deviate from 
the experiments in order for the relation (37) to be valid. 

It is interesting to notice that the change in the IR boundary conditions for the gauge 
fields with respect to the original model of [4] has some relevant consequences on the pre- 
dictions of the theory. In the original set-up the mass of the 7r(1300) resonance could not be 
reproduced and the first resonance of the axial gauge field was identified with the 7r(1800) 
state. In the present model, on the contrary, the 7r(1300) meson can be naturally accomo- 
dated. 

To check the stability of our predictions, we can compare the previous results with the 
ones obtained by using an alternative fitting procedure. For this purpose, we chose to 
minimize the largest deviation of our predictions from the experiments. In this way we 
obtained the list of result shown in the right panel of table 1. The deviations from the 
experiments are now more uniformly spread among the various observables, with a maximal 
deviation of ~ 13%. The overall RMSE is 9.6%, which is still reasonable and only slightly 
higher than the one found in the previous fit. The corresponding values of the microscopic 
parameters are L~ l = 350 MeV, M b L = 0.014, L 2 M 2 Bulk = -3.8, £L = 1.5 and M q = 
40 MeV, which are in good agreement with the previous determination. 

The heavier resonances, which have not been included in the fits, show larger deviations 
from the experimental values. For example the predicted mass for the 7r(1800) is 2140 MeV 
with an 18% deviation and for the ao(1450) it is 2070 MeV with a 40% deviation. We remark, 
however, that the heavy resonances, being close to the cut-off of the theory, are expected to 
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Experiment 


AdS 5 




Experiment 


AdS 5 




0.4 ±0.3 


0.47 


u 


0.4 ±0.3 


0.43 


L 2 


1.4 ±0.3 


0.95 




1.4 ±0.3 


0.87 


L 3 


-3.5 ±1.1 


-2.8 


L 3 


-3.5 ±1.1 


-2.4 


U 


-0.3 ±0.5 


0.0 




-0.3 ±0.5 


0.0 


L 5 


1.4 ±0.5 


0.72 


U 


1.4 ±0.5 


0.68 


Le 


-0.2 ±0.3 


0.0 


U 


-0.2 ±0.3 


0.0 


L 8 


0.9 ±0.3 


0.45 


L 8 


0.9 ±0.3 


0.39 


U 


6.9 ±0.7 


6.0 


U 


6.9 ±0.7 


5.6 


L\o 


-5.5 ±0.7 


-6.0 




-5.5 ±0.7 


-5.6 



Table 2: Predictions for the coefficients of the C(p 4 ) terms in the Lagrangian compared with 

the experimental values [18]. The values are given in units of 10 -3 . The microscopic parameters are 
fixed by the fits on the observables in table 1 (RMSE fit for the left table and 'maximal deviation' 
fit for the right table). 

have larger theoretical uncertainties. 

Other quantities that we can extract from the model are the coefficients of the 0(p 4 ) 
terms in the %PT Lagrangian, which describe the interactions of the pions with the left and 
right sources and with the spurion field related to the quark masses. The computation can 
be performed by following the holographic procedure outlined in appendix A. Due to the 
non negligible experimental uncertainty, we decided not to include these observables in the 
fit for the microscopic parameters. We also excluded from the computation the L 7 coefficient 
which arises from integrating out the Goldstone boson singlet related to the U(1)a anomaly, 
whose complete treatment is not included in the present model (see footnote 3). In table 2 
we listed the predictions of the model for two sets of microscopic parameters found with the 
RMSE and the 'maximal deviation' fit. The numerical results in the two cases are quite 
similar and show a good agreement with the experimental data. The reduced \ 2 f° r the 
RMSE fit is 1.0, while for the 'maximal deviation' fit it is 1.2, and the deviations from the 
experimental data are always below ~ 1.5 a. 

In the chiral-symmetric models without a bulk scalar field some phenomenologically suc- 
cessful relations were found among the coefficients of the xPT Lagrangian 6 : 

L 2 = 2Li , Lg = —Lio , L4 = Lq = , L3 = — 6L1 . (38) 

In our set-up all these relations remain valid, except for the last one L3 = — 6L±, which 
receives some corrections but is still well satisfied (compare [3,4]). Notice that the first and 
third relations in eq. (38), which are implied by the large N c limit of QCD [19], are not 
modified in our model. 

6 See for example [5]. 
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3 Baryons From 5d Skyrmions 



3.1 The Static Solution 

In the present model, baryons arise as 5d Skyrmions and studying their properties requires a 
slight modification of the methods of [6-8] , where the massless case has been considered. As 
a first step we will consider the static soliton configuration. A convenient and automatically 
consistent 2d ansatz is obtained, as in the massless case, by imposing the solution to be 
invariant under a certain set of symmetry transformations. These are cylindrical symmetry 
(i. e., the simultaneous action of 3-space and SU (2)y rotations), 3d parity and time-inversion, 
defined as a change of sign of all the temporal components combined with L — > —L and 
R — > —R. This leads to the following ansatz for the gauge fields 

i^(x, z) = A ± (r, z)x a x 3 + l -e ajk x k - t@. £ i*,v) &{vUi , 
< R a 5 (x,z)=A 2 (r,z)x a , ' ' (39) 

8o(x,*) = ^, 
k r 

where r 2 = ^2 i x t x % , x 1 = x l /r, e^ x ' y ^ is the antisymmetric tensor with e^ 1 ' 2 ) = 1 and the 
"doublet" tensors A^ 1 ' 2 ) are 

(40) 

Because of parity, eq. (39) also determines the ansatz for the L fields which are given by 
Z/j(x, z) = —Ri(—x.,z), L 5)0 (x, = i?5,o(— x, z) and analogously for L, R. 

The ansatz for $ is obviously obtained by imposing the same symmetries: cylindrical 
symmetry implies 

— 1 a a 

$(x,z) = u (r,z)- + iu a {r,z)— , (41) 

where 

u (r,z) = \ 2 (r,z) , 

(42) 

u a (r,z) = X 1 (r,z)x a . 

It is easy to check that parity acts on $ as $(x, z,t) — > $ t (— x, z, t), while under time- 
inversion we have $(x, z,t) —¥ cr 2 $*(x, z, —t)a 2 . Imposing $ to be invariant under these 
transformations simply implies that Ai ; 2 are real. It is useful to note that our ansatz pre- 
serves, again as in the massless case, a local U(l) subgroup of the original 5d chiral group 
corresponding to gauge transformations of the form Qr = g and g L = with 

g = exp[ia(r, z)x a a a /(2r)} . (43) 



ab 



^abc-^c 

x a x b - 5 ab 
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Under this residual U(l) the fields <p = (pi + i s and Ap in eq. (39) are respectively 
one charged and one neutral scalar and a gauge field. It is easy to check that the field 
A = Ai + i A 2 in eq. (41) also transforms as a charge-one scalar; it will be convenient to define 
its 2d covariant derivative as 

Dp\ = - tA p \ . (44) 

It is straightforward to plug the ansatz in the 5d lagrangian and to obtain the energy of 
the Skyrmion. From the gauge part in eq. (8) and (9) we obtain, as in [7,8], 



poo pz m r 

E G = 8vrM 5 / dr dz I a(z) |D f 

JO Jz vv I L 

-iL-id* \d fi {-i4?D p <t> + h.c.) + Ap P 



2\2 



1 



2 r 



(45) 



where 
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16tt 2 M 5 L ' 

while the new contribution coming from the scalar part in eq. (11) is 

„2 



(46) 



E<5> = 8ttM 5 dr dz{ a 6 (z) 



^Arp,A)-^(0A*-A0*) 2 



+ a\z)-M> Bulk X>\ 



(47) 

Notice that the total energy E = Eq + E§ does not yet give the Skyrmion mass because the 
infinite energy of the vacuum needs to be subtracted in order to get an observable quantity. 
This zero-point energy is obtained from eq. (47) by plugging in the vacuum field configuration 
which is given by 

\ = \ v = 2iv(z), (f) = (f) V =-i, (48) 

all other fields vanishing. 

The 2d EOM for the Skyrmion are easily derived, at this point, by varying the energy 
in eq. (45,47), but in order for them to be solved suitable boundary conditions need to be 
specified at the four boundaries (z = z m , z = z vv , r = and r — > oo) of our 2d space. 
At z = z vv and z = z m the boundary conditions are given, up to the sign ambiguity in 
eq. (15) that we will now fix, by the ones discussed in the previous section. The presence 
of the r = boundary merely results from a choice of coordinates, the physical 5d space 
being completely regular at r = 0. The boundary conditions will therefore result from just 
imposing regularity of the 5d fields, with no need for extra assumptions. At r — > oo we must 
require that the solution will have a finite mass, and this is ensured by imposing it to reduce, 
up to a symmetry transformation, to the vacuum configuration in eq. (48). We also want a 
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B — 1 solution, where the baryon charge B is defined in eq. (1) and is given by 



271 JO Jz m , 



B= — I <Ir I </;r'' ; ' 



dpi-ifDvt + h.c.) + Fp_ p 



(49) 



in terms of the 2d fields. The above equation can be easily rewritten (as it must, being 
the topological charge) as a Id integral on the boundaries of the 2d space, and in order to 
get B = 1 from the r — > oo boundary we must have, as in the massless case, the following 
boundary conditions 

M = (z -z ) ' (50) 

{-^m ^uv J 

s = 

which are obtained from the vacuum (48) by means of a residual U (1) transformation in the 
form of eq. (43), with a = ir(z — z vv )/(z m — z lJV ). 

Consistently, the boundary conditions for A are obtained in the same way and read 

r^oo : \ = i e Mz-z m -)/(z m -z vv ) 2v ( z y ( 51 ) 

At z — z m , the above equation implies A(r = 0, z — z m ) = —2 it;, because the "twist" 
e nr(z-z uv )/(z m -z vv ) rec | uces to — 1 at the IR while the vacuum respects the boundary condition 
(15) with the plus sign. This resolves the ambiguity and enforces the 5d Skyrmion to live in 
the minus-sign sector, with IR boundary conditions given by 

z = z m : A = -2 if . (52) 

We stress, as mentioned in the discussion below eq. (15), that the sign ambiguity in the IR 
boundary conditions results from our choice of giving generalized Dirichlet conditions on <3>, 
instead of treating it as a Neumann field and making its boundary conditions originate from 
a localized potential that would cost us more new parameters. If we had made the other 
choice, we would have had no ambiguity, and consequently no separated sectors in the field 
space. If studying this different setup in the limit of infinite strength for the coupling in the 
localized potential we expect that, while the vacuum and the meson's wave function will be 
found to fulfill eq. (15) with the plus sign, the other boundary condition will be enforced on 
the Skyrmion solution and the results of the present paper will be recovered. Coming back 
to the boundary conditions, we must still specify the ones at r = and at z — z vv . These 
are 



Ai = . _ . (z 



r^O : Q z = Zm : A = 2i ^ J M, (53) 

where the ones for r = arise, respectively, from asking the 5d field $ and its 3-space 
derivative to be regular and single-valued. For all the other fields the boundary conditions 
are the same of the massless case and are reported in appendix B. 
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3.2 Zero-Mode Fluctuations 



In order to describe the baryons we need to consider the time-dependent deformations of 
the static soliton solution. The analysis proceeds exactly as in the massless case, we will 
therefore skip most of the details and refer the interested reader to ref. [8]. 

The single-baryon states can be identified with the zero-mode fluctuations, thus an anal- 
ysis of the infinitesimal deformations will be sufficient for our purposes. The relevant con- 
figurations are the ones which describe a slowly-rotating solution, whose degrees of freedom 
can be parametrized by three collective coordinates encoded in an SU(2) matrix U(t). 

To describe the solution we need to generalize the ansatze given in eqs. (39) and (41). 
The ansatz for the gauge fields is analogous to the one for the massless pion case: 

i? A (x, z; U) = UR^x, z) , i? (x, z; U) = z) , (54) 

and 

i? (x, z; U) = UR (x, z; K)W + i Ud U^ , _R A (x, z; U) = z; K) , (55) 

where 

' i?o(x, z- K) = X{x) (r, z)k b A^ ab + w(r, z)(k ■ x)T , 

< Si(x, z- K) = [k i - (k ■ x)x l ) + B 1 (r, z)(k ■ x)& + Q(r, z)e tbc k b x c , (56) 

%(x,z;K)=B 2 (r,z)(k-x). 
The ansatz for the scalar field $ is given by 

$(x, z; U) = C/$(x, z) C/ f , (57) 

where $ is as in eq. (41) with the new definitions 

J u (r, z) = exp [-i(k ■ x)p(r, z)\ (A 2 (r, z) + i(k ■ x)r) 2 (r, z)) , 

[ u a (r, z) = exp [—i(k ■ x)p(r, z)\ [(Ai(r, z) + i(k ■ x)rji(r : z)) x a — i ((k ■ x)x a — k a ) ((r, z)\ . 

(58) 

In the above equation, the 3-vector k a denotes the Skyrmion rotational velocity 

K = k a a a /2 = -irfdU/dt , 

and the field p is the same that appears in the ansatz for the gauge fields in eq. (56). 
The above ansatz can be obtained, similarly to the one for the static solution, by imposing 
time-inversion, parity, and cylindrical symmetry. 
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Plugging the ansatz in the 5d lagrangian one obtains the collective coordinates lagrangian 



L = —M + ^ k a k a , 



(59) 



where M is the Skyrmion mass and A is its moment of inertia. The latter receives a contri- 
bution from the gauge Lagrangian 

X G = 16nM 5 ^ / dr / dz\a(z) 



■y poo rziR 
'3 / dr 

° JO Jz m 



[Dppf - r 2 (dpQ) 2 - 2Q 2 - 



+r 2 {D n xf + y {dfiw) 2 + (x{x)X{x) + w 2 ) (l + ^(^^(a,)) - 4w<f) {x) X(, 

+1 L - toPDppxto (Ax% } + 2c^9 p (rQ) x^)^ (D p( f>) {y) 

-w Qe^Spp (^ (ie) - 1) + r Qe^Ap^j + 2r Qe^Dppd, (^) | 

and a contribution from the scalar, which is given by 
16 



(60) 



A 



= — 7rM 5 / dr dz \ a 3 ( :) 
3 Jo iz uv 



jX'XiDppXDpp) - ^ (Dpi])* (Dpi]) - y(fyCXfyC) 



[A(£>^)* - ^A)* + h-c] - \v*V + \W ~ H? ~ j(x^ ~ Ax*) 



*\2 



-<(#*-0?f)--C 2 (i + 0*0) 



r 2 a 5 (z)M 2 Bulk 



1 * 1/2 



where we defined 



ri = i] X + iri2, Dpi] = dpi] - iAprj . 

while the other notations are defined in appendix B. 



(61) 



(62) 



3.3 Numerical Results 



The soliton solution can not be determined analytically, however, it can be studied numer- 
ically by using the techniques described in [8]. In the massless case it was found that, due 
to the peculiarity of the 5d gauge action, the soliton solution is stabilized thanks to the 
presence of the CS term [7]. This peculiar feature disappears once we modify the action 
by the introduction of the bulk scalar field and, in the present model, we checked in our 
numerical analysis that the Skyrmion size is stable even if the CS term is not present. 

From the soliton solution we can extract the electromagnetic and axial properties of the 
nucleons, which are encoded in a set of form factors which parametrize the matrix element 
of the currents on two nucleon states. 
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The chiral currents can be determined by computing the variation of the action with 
respect to the sources 1^ and r M . It is simple to show that the action describing the scalar 
field $ does not contribute to the currents, which are exactly the same as in the massless 
pion case: 

J a Lti = M 5 (a(z)L^ 5 ) \ z=Zm , J L , = M 5 (a{z% 5 ) \ z=Zvv , (63) 
and analogously for R. 

The isoscalar and isovector form factors are defined through the relations 

(N f ffl2)\J°s(0)\Ni(-?m) = G s E {f) X \xi, 

(N f (q/2)\4(0)\Ni(-q/2)) = < ^jj^X/ 2 ^ X » 
(N f (q/2)\4 a mNi(-m) = G V E (q 2 )x\(2I a )Xi, 

(N f m)\Jy(0)m-q/2)) = < X & ( 2/ °) ^ > ( 64 ) 

where the currents are given by J v = + J£ and J s = 1/3(Jr + Jl), and we used the 
notation Xi,f f° r the nucleon spin/isospin vectors of state (normalized to x^X — 1) an d the 
definition (S x q) 1 = e tjk S j q k . From the axial current J% = Jr~ Jb we define the axial form 
factors 

(N f m)\JT(0m(-q/2)) = x\[^G A (q 2 )S T + (G A (q 2 ) - J^tf 2 )^l]/^ , (65) 

(7V /(g7 2)| J}«(0)|^(-9y2)> = (66) 

where St = S — q S • q and Sl = q S ■ q are the transverse and the longitudinal components 
of the spin operator. 

To find the explicit expressions for the form factors we need to plug the ansatze for the 
soliton solution into the definitions of the currents (63) and then perform the quantization 
of the soliton solution as explained in [8]. The result is the same as in the massless pion 
case: 

N f 

G S E = --^l J drrjo(qr) (a(z)d z s) uv 

G e = ^J drr 2 Jo (qr) [a(z) (d z w - 2 (D zX ) {2) )~ 

<% = ?^J*r^(a(z)d.Q) uv 
v M N N C f 2 ji(qr) / \ 

M = J drr -jr \< z ) ( D ^h)) 

G A = J drr a W~~ ~ rA z^j - a(z)(D z (f)) {1) jo(qr) 



uv 



uv 

(67) 

uv 
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Experiment AdSs Deviat. Experiment AdSs Deviat. 



M N 


940 MeV 


1090 MeV 


16% 




0.44 


0.43 


2% 


II T r 


2 35 


1 18 


100% 

J- \J\J / (J 


9a 


1.25 


0.58 


100% 




0.79 fm 


0.82 fm 


4% 




0.93 fm 


0.97 fm 


4% 


\J ( T M,s) 


0.82 fm 


0.84 fm 


2% 


\J ( r M,v) 


0.87 fm 


0.87 fm 


0.5% 




0.68 fm 


0.65 fm 


5% 



M N 


940 MeV 


1104 MeV 


17% 


Us 


0.44 


0.43 


3% 








i nn% 


9 A 


1.25 


0.59 


100% 




0.79 fm 


0.84 fm 


6% 




0.93 fm 


1.02 fm 


9% 


\J ( r M,s) 


0.82 fm 


0.86 fm 


5% 


\j ( r M,v) 


0.87 fm 


0.86 fm 


1% 




0.68 fm 


0.68 fm 


0.2% 



Table 3: Prediction for the static nucleon observables with the parameter values fixed by the fit 
on the mesonic observables (RMSE fit for the left table and 'maximal deviation' fit for the right 
table) . 



where j n are spherical Bessel functions. 

By employing suitable numerical techniques, the 2d EOM 7 obtained by varying the 
soliton mass M and its moment of inertia A can be solved, and both the static and slowly- 
rotating Skyrmion solution computed. The numerical predictions for the static nucleon 
observables are listed in table 3. In the analysis we used the values of the microscopic 
parameters obtained from the fits on the mesonic observables presented in section 2.2. 

The numerical results for the two sets of microscopic parameters considered show very 
similar deviation patterns from the data. Many of the numerical predictions are very close 
to the experimental results, although the magnetic vector moment fiy and the axial coupling 
qa present a deviation of order 100%. We notice, however, that the overall agreement with 
the data (root mean square error 45%) is still compatible with the possibility of having 
sizable 1/N C corrections, which can not be excluded given that N c = 3. 8 By comparing the 
present results with the ones found in the simplified model without a pion mass [8], we see 
that all the observables show an improved agreement with the data except for the fiy and 
qa, whose deviations become significantly larger. 

Form the qualitative point of view, we remark that an important check of the validity of 
7 The EOM for the 2d fields are reported in appendix B. 

8 It is interesting to notice that using a different approach to the quantization of the collective coordinates, 
as suggested in [20], one gets much better predictions for /iy and gA- With this procedure, the predictions 
for [i v and gA are rescaled by a factor 5/3, thus agreeing with the data at the 20% level, while all the other 
observables are unchanged. We have no reason to believe that the modified quantization procedure correctly 
captures the 1/N C corrections, nevertheless, this result seems to point out that large corrections could indeed 
be responsible for the deviations of \iy and gA- 
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the description of baryons as solitons is the behavior of the electric and magnetic vector radii, 
namely te,v and tm,v- These two quantities are expected to be divergent in the chiral limit, 
as explicitly verified in [8], while they should become finite once a pion mass is introduced, 
as we find in the present model. 

4 Conclusions 

We have shown that it is rather easy to construct a model of holographic QCD which 
describes at the same time the pion mass, the QCD anomalies and the baryons as topological 
solitons. After introducing an explicit minimal model we have studied its phenomenology in 
both the mesonic and baryonic sector and found a significant level of agreement. In extreme 
synthesis, our result is that the general picture on the holographic QCD models outlined in 
the Introduction survives unchanged to the inclusion of the pion mass. 

Few unexpected results have been found, however, that is worth discussing. In Sect. 2 
we saw that our model easily reproduces the mass of the 7r(1300) meson, a task that was 
impossible to achieve in the original scenario [4]. This came because of the change in the 
IR boundary conditions and gives a phenomenological support to this modification, that 
we had instead motivated on purely theoretical grounds. It is also remarkable that the 
other predictions are almost unaffected so that all the valid phenomenology of the original 
construction is retained. The "new" observables that were absent in the original model, i.e. 
the anomalous parity couplings originating from the CS term, also show a good agreement 
with the observations. Our results in the baryonic sector, shown in table 3, are also surprising, 
especially if compared with the ones obtained in the chiral-symmetric case [6-8]. For all the 
observables except \iy and qa, a significant improvement is found in the agreement with 
observations. The isovector radii, that have became finite due to the presence of the pion 
mass, are also extremely well predicted. The situation has got significantly worst, on the 
contrary, for \iy and g& that have became a factor of 2 smaller than the observations. 9 This 
failure persists in both the best fit points that we have used as input parameters in table 3, 
so that it is probably a robust feature. It might signal that the model is incomplete, but it 
might also be attributed to anomalously large 1/N C corrections. 

Some final comments on the theoretical implications of our results. The microscopic 
origin of holographic QCD models is basically unknown, even though the holographic imple- 
mentation of the chiral symmetry provides a robust (but purely technical) connection with 

9 Actually, the discrepancy is almost exactly given by a factor of 2. By looking at table 3 one could 
imagine a factor 2 mistake in the problematic predictions, we are however confident of our calculations. 
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AdS/CFT. The success of the Regge phenomenology suggests, independently of AdS/CFT, 
the dual of large-iV c QCD being a string model and the validity of the holographic QCD 
approach suggests that this string model should contain a sector that is well described by a 
5d field theory similar to the one we have considered. In the case of exact chiral symmetry, 
the Sakai-Sugimoto model [21,22] provides a partial realization of this idea because it is 
equivalent to holographic QCD for what the physics of the light meson is concerned. 10 It 
is on the contrary different, and problematic, in the baryonic sector [6-8, 12] and its phe- 
nomenology in the sector of higher spin states (where it genuinely differs from holographic 
QCD and shows its stringy nature) seems not very promising [22]. The inclusion of the 
explicit chiral breaking considered in the present paper provides an additional piece of infor- 
mation. The chiral symmetry is not explicitly broken in Sakai-Sugimoto, and even though it 
was possible to include the explicit breaking by some deformation, the resulting model could 
not reduce to a field-theoretical model such as the one we have considered. The reason for 
this is that the Left and Right global groups are localized, in Sakai-Sugimoto, at two dif- 
ferent boundaries of the 5d space and the quark mass spurion M is unavoidably a non local 
object which is impossible to describe in a field theoretical language. Therefore, any string 
model that incorporated ours, inheriting its phenomenology, would be deeply different from 
Sakai-Sugimoto; it might be worth trying to construct one following a bottom-up approach. 
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A The Effective Action for the Pion 

In this appendix we will compute the effective action for the pion at 0(p A ). Given that a 
complete treatment of the U(1)a anomaly is not included in the model, in the following we 
will neglect the Goldstone boson related to this symmetry, namely the rj' field, and we will 
only consider a model with a chiral symmetry SU(2) L x SU(2) R . 

An efficient way to perform the computation is to use the holographic approach presented 
in [14]. At tree-level, the holographic action for the pion is given by the 5d action for the 

10 Scc [10] for a precise justification of this equivalence. 
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gauge and the scalar fields (eqs. (8) and (11)), where the 5d fields satisfy the bulk EOM's 
with the usual IR boundary conditions given in eqs. (14) and (15). The UV conditions are 
modified as 

L,\ Zvv = U(x)(l, + id,)U(xy 

^Lv =r V A , (68) 

*L = tf(*)(£) 

where U (x) represents the 4d Goldstone matrix, which transforms as 

U(x) g R U{x)g{ (69) 

under a chiral SU(2)l x SU(2)r 4d transformation. Notice that we are not interested in 
possible terms involving the scalar and pseudoscalar sources, so we did not include any source 
term in the UV condition for the scalar field in eq. (68). 

To derive the complete effective action for the pion one would need to solve the full 
bulk EOM for the 5d fields. However, due to the presence of interaction terms, this can be 
done only perturbatively. As usually done in xPT, we use an expansion in powers of the 
momentum p and we treat the external sources and r M as 0{p) terms, while M q will be 
treated as an 0(p 2 ) term. 

We expand the solutions of the EOM using a mixed momentum-space representation 
V, = ft{z)%{p) + V^\p, z, V, A, U, Mg) , 

A, = f° A (z)Mp) + A®(p, z, V, A, U, M q ) , (70) 
$ = $(o)( z ) + $(2)( p , z , v, A, U, M q ) + $( 4 )(p, z, V, A, U, M q ) , 

where we denoted by V and A the values of the vector and axial gauge fields at the UV 
boundary. Notice that, due to the tensorial structure, the gauge field solutions start at 0(p) 
and their next to leading terms are of 0(p 3 ), while the scalar field can be expanded in terms 
of 0(p 2n ). All the terms in the expansion of the gauge fields satisfy the same IR boundary 
conditions as the original fields: 

d z V ll {x,z m ) = Q, ^ 
Au(x,z m ) = 0. 

Instead, for the scalar field we have 

f *<»)(*,,„) -€1. 

[$^{x,z m ) = fori>2. 

The UV boundary conditions for the gauge fields are chosen so that the higher terms in the 
expansion vanish, while for the leading terms we have 
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For the scalar field we impose 

f $(°)( 2W ) = 0, 
<£ (2) (^v) = t/(^) ~M 9 1, (74) 
$( 4 )( Zuv ) = fori>4. 

By using the bulk EOM's and the boundary conditions for the fields, one can verify that 
the terms of 0(p 3 ) in the expansion for the gauge fields and the ones of order 0{p 4 ) in the 
expansion for the scalar field do not contribute to the effective action for the pion at 0(p A ) 
(see the discussion in [14]). 

As a first step of the derivation of the effective action, we will consider the contributions 
coming from the 5d gauge action in eq. (8). It is convenient to rewrite eq. (8) in the form 

S g = -2M 5 f d 4 x a(z lJV ) Tr [V^d z V + A^d z A»] 
Juv 

-M 5 J d 4 x J dza(z) Tr [L^L^ + R^R^] . (75) 

The leading terms in the expansion of the 5d fields satisfy the equations 

d z (a(z)d z f°(z)) = 0, 

d z (a(z)d z f A (z)) - 2a\z)^ff A (z) = , (76) 
8 Z (a 3 (z)d z &°\z)) - a 5 (z)M*$W)(z) = 0. 

The solution for the $^°- ) field can be obtained from eqs. (17) and (18) by setting M q = 0. 
The equation for the vector gauge field admits the simple solution fy(z) = 1, while the 
equation for f\ in general can not be solved analytically. 

By substituting the above expressions in the gauge action and using the relation 



A, = l -U (D.l/y = -- (D,U) £/+ = u M , (77) 

we find that the first line in eq. (75) gives the kinetic term for the pion, from which we can 
extract the pion decay constant 

fi= -2M 5 a(z)dJ° A (z)\ z=Zvv . (78) 

From the second line in eq. (75) we get contributions to the 0(p 4 ) terms in the pion effective 
Lagrangian. In the standard form of %PT [19] we get the following contributions 

L?> = ^ [dza(z)(l-(f A rr \f = 2L ?, 

M M\ f { 4 9) = -QL[ 9) ■ (79) 
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Now we consider the contributions to the effective action coming from the scalar action 
in eq. (11). To derive the action we need to compute the term in the expansion of the 
5d scalar field. This term satisfies the bulk EOM 

±-8 z (a\z)d z & 2) ) + a\z)Ml^ = -2$(°> [-2(/»)V - </S^«1 » ( 80 ) 



cfi(z) 

where we defined = <9 M w M — i\V^ u^]. The solution can be split into two parts: 

* (2) = *ff + *?> , (81) 
where $^ is a solution of the homogeneous part of eq. (80) with boundary conditions 

= ^ (^) ^ M 9 1 , §2W = , (82) 

while ^q 2 - 1 satisfies eq. (80) with boundary conditions 

<S>£\z lJV ) = <S>£\z m ) = . (83) 



The solution for $^ is simply given by eqs. (17) and (18) with the choice £ = 0. The 
solution for $q ^ can not be found analytically, and we will parametrize it as 

<$>T = f { £\z)u^ + f¥\z)DX- (84) 

From the action for the scalar field we get an 0(p 2 ) contribution to the pion effective 
action, which can be written as 

=M 5 f d*xa 3 (z)Tr \<f>$ d z & 0) + h.c.l . (85) 
Juv L J 

To obtain this expression we integrated by parts the terms containing derivatives with respect 
to z and we used the bulk EOM and the boundary conditions for the terms in the scalar field 
expansion. Eq. (85) corresponds to a mass term for the pion, which in the limit z m — > 
becomes 

= 2^ a £ J d 4 x Tr [UM\ + M q tf] . (86) 

The expression for the pion mass and for the pion decay constant can be easily computed as 
an expansion in the parameter £. Some approximate expressions are reported in eq. (36). 
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Computing the 0{p A ) terms in the action we get 

Sf ] = M 5 J d 4 x J dza 3 (z)[2Tr [$(°> ff\z)f A {z)D^D v u v 

+4Tr [$(°)/« {z){f A {z)) 2 uX^u 
+iTr ' 
+2Tr 



* (0) ($2 }t - *2 } ) /M^x] 



+^ /" d 4 xa 3 (z)ivk(/f (^uX-i/fW^XJ^+h-c-l • ( 87 ) 
^ Juv L J 

The 0(p 4 ) action for the pion can be simplified by using the EOM for the pion field coming 
from the 0(p 2 ) effective action. Using the standard notation of xPT [19], the kinetic and 
mass terms for the pion are written as 



C = & {Tr [(D.UyD^U] + Tr ( X *U + X U ] ) } ■ 
From this Lagrangian we get the EOM for U : 

(D^U)^ - U(D^Uy + U X ] ~ = , 



(88) 



(89) 



which can be rewritten as 

AiD^ = X U ] - Ux ] ■ (90) 

Moreover, by comparing eq. (88) and eq. (86), we can extract the relation between M q and 
X, which, in the z vv — > 0, limit reads 



M q = —^ X 
q M 5 a£ 8 A 



(91) 



By using these relations we get the scalar contributions to the coefficients of the 0(p 4 ) 
effective pion action 



' Lf = f \dz a 3 (,)Tr (*(°>) fl\z){f a {z)f , 



KM: 



a 3 (z)f£\z)d,fg\z) z v +2j dz a%z)f$(z)(f a (z)fTr{<S>^)j , 

^ ^(z)fS\z)djr\z)\ Zv -J dz a\z)f a (z)Ur (z) -2Kf£\z))Tr (*<°>)} , 

(92) 

where k is defined by the relation M q = k X . 
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B The Equations of Motion 



In this appendix we report the EOM for the 2d fields which appear in the ansatz for the 
zero-mode soliton fluctuations in eqs. (39), (56) and (58), and we summarize the notation 
used in the paper. 



B.l The Equations of Motion 

Before writing the EOM's for the 2d fields, it is useful to recall the residual symmetries 
which survive after we choose the ansatze for the Skyrmion solution. As already discussed 
in sec. 3.1, the ansatz preserves a U(l) local symmetry with 2d gauge field A^, which cor- 
responds to the 5d gauge transformations given in eq. (43). The fields 0, A, x an d i] are 
charged under this symmetry, thus it is convenient to write the action and the EOM in terms 
of their covariant derivatives 

D^X = <9 M A - iA^X , D^r] = d^r] - iA^r] . 

There is also a second residual U (1) generated by the U (1) l,r 5d transformations of the form 
gn = g and g L = g^ with 

g = exv[0(r,z)(k-x)\ , (94) 

whose associated gauge boson is B p (see [8]). The p field transforms as a Goldstone boson 
under this symmetry, so we can define its covariant derivative as 

D^^d^p-B^. (95) 

The EOM's for the 2d fields can be easily found by substituting the ansatz for the 
Skyrmion solution into the 5d action. The EOM's for the fields which appear in the static 
soliton case 

' D* (a(z)D-^) + ^0(1 - |0| 2 ) + ^A(A0t _ A t) + ryLe^fy (£) D D <P = 

& {r 2 a{z)A flD ) - a(z) {i^D 9 <t> + h.c.) - l -a\z)r 2 [X^(DpX) - (^A)U] 

+ 7 ^%(f)(|0| 2 -l) = O > ( 96 ) 

d fi (a(z)d^s) - [d^-itfDvt + h.c.) + A^} = 

k (r 2 a 3 {z)D n X) - a 3 (z)(f)(Xft - 0A+) - a 5 {z)r 2 M 2 Bulk X = 
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while the equations for the fields which are turned on for the rotating Skyrmion are 

1 



d^{r 2 a{z)df i w) - 2a{z) [w(l + |0| 2 ) - xtf - </>X f ] + iLe^ 



1 



D\r 2 a{z)D vX ) + a(z) [2w<P - (1 + |0| 2 ) X ] + -a\z)r 2 \(\ X ] - X^) 



- -yLe&iDtf) \id D {rQ) + D p p] = 
-d~^r 2 a(z)d fl Q) - -a(z)Q 



4 LV 



{iD-^{D, X ) ] + h.c.) + l -A-^{2w - X <P j 



a 



2 n 

Dfipdp 



d fl (a(z)D fl p) - -a 3 (z) [(tjA* + - 2<(0A f - A^)] 



— e 



{D-^{D dX ) ] + h.c.) + -A-^x ] - xtf) 



a 



2 d p (rQ)d 9 ( 



d p {r 2 a{z)B Dp ) + 2a(z)D ft p + -r 2 a 3 (z) [2\^\D p p + {r]{D-^ - \{D-^ + h.c.)] 

+ 7 Le^{ [(x - wct>){D p ^ + h.c] + (1 - |0| 2 )«9 P «; - d p (^) } = 

D^(a 3 (z)r 2 D f]i r]) - \d^(a 3 (z)r 2 D n p) - 2a 3 (z)r 2 (D»p)(D fl \) 

- a 3 (^(2 + |0| 2 ) + a 3 (z)0(077t + 4<) - r 2 a\z)M 2 Bulk r ] = 

[ ^(r 2 a 3 (z)^C) - a 3 W [C(l + I0| 2 ) " *(<^ f - # f )] - r 2 a 5 (^)Ml uifc C = . 



(97) 



In order to find suitable equations for the numerical analysis of the solutions, the EOM's 
must be rewritten as a system of elliptic partial differential equations. For this purpose 
we need to choose a gauge fixing condition for the residual 2d U(l) gauge symmetries. A 
possible choice is a Lorentz gauge condition 



9M, = 0, 



9^ = 0. 



(98) 



With this condition the equations for A v become J v = dp (r 2 aA> IU ) = r 2 adp i d pl A v +dp i {j- 2 a)A^ v 
which is an elliptic equation and a similar result is obtained for B^. 



B.2 The Boundary Conditions 

The derivation of the boundary conditions for the 2d fields has been discussed in section 3.1. 
Here we report the list of conditions we need to impose on the scalar field components as 
well as the conditions for the gauge fields, which are analogous to the ones for the massless 
case [8]. 

The IR and UV boundary conditions on the 2d fields follow from the boundary conditions 
for the 5d fields (eqs. (14) and (10) with vanishing sources for the gauge fields, eqs. (15) and 
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(6) for the scalar) and from the gauge choice in eq. (98). They are given explicitly by 



z = z n 



and 



Z = Z v 



0i = O 

d 2 <t>2 = 
Ai = 
d 2 A 2 = 
d 2 s = 
I A = -2i£ 



Ai = 
<9 2 A 2 = 
s = 

k A = 2iM g (z uv /2: IR ) A ~ 



Xi = 
<9 2 X2 = 
d 2 w = 
d 2 Q = 




p = 
Bi = 
<9 2 £ 2 = , 
?7 = 

C = o 



(p = 

B 1 = 
d 2 B 2 = 
T] = 

I C = o 



(99) 



(100) 



The boundary conditions at r = oo are obtained by the requirement that the energy 
of the solution minus the vacuum energy for the scalar field is finite. To obtain a soliton 
solution with B = 1 one imposes 



oo 



( _ _j e in(z-z vv )/(z m -z m ,) 

<9iAi = 

A 2 = ii I (z m - z vv ) 
s = 



■y^ — j^^ n ( z ~ 2: uv)/(^ : m — Z\-v) 
W — — 1 

Q = 



fp = 

= 
B 2 = . 

T] = 

C = o 

(101) 



On the r = boundary of the domain we must require the 2d solution to give rise to 
regular 5d fields and the gauge choice in eq. (98) to be fulfilled. These conditions are fulfilled 
with the choice 



r = : 



(l + 2 )/r ^0 
A 2 = 
9iAi = 
s = 
Ai = 
diA 2 = 




' p/r ->■ Bi 
= 
B 2 = 

ifc = C 

772 = 

I <9iC = 



(102) 



C The QCD Anomaly 

In this appendix we describe different forms of the QCD anomaly and discuss the relation 
with the CS term included in the 5d theory. Part of the material that follows overlaps with 
appendix A of [10]. 
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The CS term (9) can be written as 



Scs — 



24vr 2 



J d 5 x[uj 5 (L)-UJ 5 (R)} , 



(103) 



where U5 differs by a total differential from the standard text-book CS form: 



w 5 (A) = % Tr 



A (dA) 2 + 3 A 3 dA + 3 A 5 
2 5 



= 1 ^ A(dA) 2 + i - ATi [F 2 ] + 1 d 



ATr 



AF A 3 

4 



cu 5 (A) + dX(A) . (104) 



The variation of the CS is given by eq. (13), where the 4-form 

u>\(a,A) = X -a [dA) + ^STr [F 2 ] , (105) 
is defined from the relation 5 a oJ5 = duj\ and it is related to the standard u\ by 

uj\(ct,A) =u\(a,A) -5 a X(A), (106) 



where 



ul(a,A) = Tr 



ad ^AdA + ^A 3 



(107) 



Provided the IR term in eq. (13) is cancelled, the CS variation gives the anomaly of 
eq. (4), which however does not coincide with the standard text-book QCD anomaly that is 
normally put in the "symmetric" form 



N c f 



(108) 



The two forms of the anomaly (A and A sym ) are equivalent because they only differ by a 
local counterterm 



A = A - 



[6 aL X(L) - 5 aR X(R)] 



(109) 



sym 247T 2 

Obviously we may equivalently have added the X local counterterm to the 5d Lagrangian 
and kept the standard form of the QCD anomaly. This would have not affected any of our 
results because X only depends on the 4d 1 and r sources, whose physical value is zero. 

Similarly, the QCD anomaly could be also put in the Adler-Bardeen form. Starting from 
the symmetric anomaly, this is achieved by the addition of the Adler-Bardeen counterterm, 
as explained in [14]. 
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